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Abstract 

In this paper, we formulate the relativistic heat equation and the relativistic 
kinetic Fokker-Planck equations into the GENERIC (General Equation for 
Non-Equilibrium Reversible-Irreversible Coupling) framework. We also show 
that the relativistic Maxwellian distribution is the stationary solution of the 
latter. The GENERIC formulation provides an alternative justification that 
the two equations are meaningful relativistic generalizations of their non- 
relativistic counterparts. 
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1. Introduction 

In this paper, we formulate the relativistic heat equation and the relativis¬ 
tic kinetic Fokker-Planck equations into the GENERIC (General Equation 
for Non-Equilibrium Reversible-Irreversible Coupling) framework and discuss 
its consequences. Let us start by recalling some properties of the classical 
heat equation 


d t p = is A p, 

and the kinetic Fokker-Planck (Kramers) equation, 


( 1 ) 
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Eq. (PP) describes the distribution of heat (or variation in temperature) in 
space over time, where v > 0 is the thermal diffusivity. In Eq. ([2]), the spatial 
domain is M? d with coordinates ( q,p ), with q and p each in W l . Subscripts 
as in divg and div p indicate that the differential operators act only on those 
variables. The function V = V(q) is given, as are the positive constants m, 7 
and 9. This equation characterises the evolution of the probability density 
of a Brownian particle moving in a homogeneous, viscous medium and under 
the influence of external force fields. In this context, q,p and m represent 
the position, momentum and the mass at rest of the particle; 7 is the friction 
coefficient and 9 = kT where k is the Boltzmann constant and T is the 
absolute temperature. 

It is well-known that both Eq. (JT]) and Eq. ((2]) allow infinite speed of 
propagation. However, this violates Einstein’s theory of special relativity 
that the speed of light c is the highest admissible velocity for transport of 
radiation in transparent media. To remedy this issue, generalizations of Eq. 
dU) and Eq. (J2J) to the relativistic setting are required. For the heat equation, 
it was carried out first by Rosenau in [ 3 ^] and later on by Brenier in [0], 
Their idea is to change the classical flux based on Fourier’ law, F = zA7p, by 
a flux that saturates as the gradient becomes unbounded, leading to models 
of nonlinear degenerate parabolic equations, and in particular the following 
relativistic heat equation 


dtp = v div 




y + 72 


|Vpp 


(3) 


where c is the speed of light. Eq. (J3]) and related models have been investi¬ 
gated further by many authors, see e.g., [i|, [ij and references therein. 



authors proposed the following equation as the relativistic generalisation of 
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d t p = - div„ 


cp 


\Jm 2 c 2 + \p\ 2 


P 


+div. 


v q v + T 


cp 


■\Jm?c 2 + \p\ 2 


P 


In 00 the authors suggested an alternative one, which is 


+^9A p p. 

(4) 


d t p = - div„ 


cp 


y/m 2 c 2 + |p| 2 


P 


+ divp 

{y q V + 7 —) P 


LV H mJ J 


7 # divp 


m c 


\Jm 2 c 2 + |p | 2 


/ + 


p <g)p 

m 2 c 2 


X7 p p 


(5) 


Notations in these equations are similar to those in Eq. (]2j). Additionally, 
/ is the d- dimensional identity matrix, p<g)p denotes the d- dimensional matrix 
with entries ( p®p)tj = , and c is the speed of light. The difference between 

Eq. (J4]) and ((5]) lies in the last two terms on the right-hand sides, namely 
the drift terms in p— directions and the diffusion matrices. There is a large 
literature on these two equations, see e.g., 0, [ll|, 13, 17, 25, 26, 27, 0, 0] or 


review papers 0,01 and references therein. Many properties of them, such 
as ih-theorem, fluctuation-dissipation relation and stationary solutions, have 
been studied in the cited papers. More particularly, in [g] and 003 
the authors have explained why (J3J) and ([5]), while having finite speed of 
propagation as desired, are meaningful relativistic generalizations of (|T|) and 
([2)) respectively. The argument is that one can recover (JT]) and ([2]) respectively 
from ([3]) and ([2]) as the speed of light tends to infinity, i.e., c —* oo. In 
addition, in 0, 0, 3, the authors suggest that Eq. (J5]) is more physically 
preferable than other models derived in 0] because it has the relativistic 
Maxwellian distribution as the unique stationary solution. 

The aim of this paper is to formulate the relativistic equations, (E]), (J1D 
and (|2]), into the GENERIC framework, which is a well-established frame¬ 
work for non-equilibrium thermodynamics. GENERIC is a formalism for 
non-equilibrium thermodynamics unifying both reversible and irreversible 
dynamics. It has been used widely in mathematical modelling of complex 
systems. It is now nearly impossible to list all references abouh GENERIC, 


so we will mention some recent ones. In the original papers [24], [38], it was 


Tn 0 , the external potential V = 0. 
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originally introduced in the context of complex fluids with applications to the 
classical hydrodynamics and to non isothermal kinetic theory of polymeric 
fluid. Recently it has been applied further to anisotropic inelastic solids 28], 
to viscoplastic solids j29|, to thermoelastic dissipative materials [32|, to the 


soft glassy rheology model [1 
Timoshenko and damped Bresse systems 
tic setting have been investigated in, e.g., 


to turbulence 37] as well as to the damped 
2ol l. GENERIC in the relativis- 
Research on GENERIC 


from rigorously mathematical perspectives is actively carried out, see e.g., 
32, MM £ • 

Placing (l3|) . (JU) and ([5]) in the GENERIC framework have a couple of 
benefits. Firstly, let us remind that equations ([I]) and ([2]) have already been 
cast into the GENERIC setting in [l| and [l8[. We show that the GENERIC 
structure of (j3j) , and of (J4J) or ([5]) are analogous to that of ([Tj) and (J2]) re¬ 
spectively, thus putting all of them in a common framework and providing 
an alternative justification that (]3]) and (j4]) or ([5]) are meaningful relativistic 
generalizations of ([T]) and ((2]) respectively. A crucial property used in the 
construction is a fluctuation-dissipation relation which holds true for both 
non-relativistic and relativistic models and is more general than the one in 
12, Eq. (13)], see Remark |T] for more discussion. We also show that, as a 
consequence of the GENERIC formulation, (J4]) and (J5]) have the same sta¬ 
tionary solution, which is the relativistic Maxwellian distribution. It should 
be mentioned that stationary_solutions of Eq. (HJ) and Eq. ([5]) have been 


studied in the literature [12|, |14|, [2[. Here we provide an alternative computa¬ 
tion using the GENERIC formulation. Our approach sheds light on this issue 
because we show that Eq. ()4]) and (J5]) actually have the same form, see (l26j) 
below, and that’s why they share the same stationary distribution. Secondly, 
as being GENERIC systems, (J3]) and (J2]) will automatically justify the first 
and the second law of thermodynamics. Thirdly, the GENERIC framework 
reveals physical and geometrical structures of ((3]) and (J5]) by specifying ex¬ 
plicitly the conservative part, the dissipative part, the energy functional and 
the entropy functional as well as their role as driving force of the two parts 
respectively. The reversible-irreversible splitting structure of the relativistic 
kinetic Fokker-Planck equations seems unclear in the existing papers. Last 
but not least, we stress that our derivation is applicable for arbitrary dimen¬ 
sion d, not necessarily to separate into 1-dimensional and 3— dimensional 
cases as in 0 0. 

The rest of the paper is structured as follows. In Section [2] we summarize 
the GENERIC formalism in general and of ([I]) and ([2D as shown in [l] and 18 
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respectively. We formulate the relativistic heat equation and the relativistic 
kinetic Fokker-Planck equations in the GENERIC framework in Section 13.11 
and in Section 13.21 respectively. In Section 13.31 we compute the stationary 
solution of ()5|) . Finally, a conclusion and discussion of further development 
is given in Section [4] 


2. GENERIC framework and its generalisation 


2.1. GENERIC framework 

In this section, we summarize the GENERIC framework. We refer to 
the original papers [24j, 38| and the book 36] for further information. To 
formulate a thermodynamic system in the GENERIC framework one needs to 
specify five building blocks {Z, L, M, E, S} and verify certain conditions put on 
them as following. Let z G Z be a set of variables which appropriately describe 
the system under consideration, and Z denotes a state space. Let E, S : Z —> M 
be two functionals, which are interpreted respectively as the total energy and 
the entropy. For each z 6 Z, let L(z) and M(z) be two operators, which are 
called Poisson operator and dissipative operator respectively, that map the 
cotangent space at z into the tangent space at z. Associated with the Poisson 
operator and two observables, F and G (i.e. real valued, sufficiently regular 
functionals of the set of variables z), is a Poisson bracket 


{F , G}L: .ffih. L(z) ^w 


8 z 

satisfying an anti-symmetric condition 


8 z 


( 6 ) 


{F.G} l = -{G,F} l , (7) 

and the Jacobi identity 

{{F i, F 2 } l , F 3 } l + {{F 2 , F 3 } l , FJ l + {{F 3 , FJl, F 2 } l = 0, (8) 

for three observables F j , F 2 , F 3 . Similarly associated to the dissipative oper¬ 
ator M = M(z) and two observables, F and G is a dissipative bracket 


|F ,G] M :.TWW G « 


8 z w 8z 

which is symmetric and positive semi-definite, i.e., 

[F. G]m = [G, F]m, [F,F]„ >0. 


(9) 


( 10 ) 
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Moreover, { L, M, E, S} are required to fulfill the degeneracy conditions: for 

all z G Z, 


LM^Uo, 


M(z 




= 0. 


( 11 ) 


<5z v ' dz 

A GENERIC equation for z is then given by the following differential equation 


da = L(z 


*E(z) 

dz 


M(z 


SS(z) 

5z 


( 12 ) 


where da is the time derivative of z; || are appropriate derivatives of E 
and S respectively. They could be either the Frechet derivative or a gradient 
with respect to some inner product. Respectively, the meaning of the dot in 
(ED and © is that of the duality pairing or the formal inner product. 

The degeneracy conditions (TTTT) together with the symmetries of the Pois¬ 
son and dissipative brackets ensure that the energy is conserved along so¬ 
lutions of (TT2|) and that the entropy is a non-decreasing function of time. 
Indeed, by straightforward calculation, we have 


and 


dE(z(t)) 5E(z) 


dt 


d z 

GD <5E(z) 


dz {nj AE(z) 

dt 


dz 


L(z 


*E(z) 

'~dT 


Mfz 


8S(z) 

dz 


dz 


L(z 


^E(z) 

'~et 


i{E,E} L i 0 , 


dS(z(t)) dS(z) 


dt 


dz 

rm <5S(z) 


dz (p d S(z) 
dt dz 
dS(z) 


L(z)^ + M( 

dz 


, *S(z) 

dz 


dz 


■ M(z 


dz 


jm m 

= [S, S] M > 0. 


Moreover, GENERIC allows us to determine the stationary solution. Equi¬ 
libria are the maximizers of the entropy S under the constraint that energy 
is constant E(z) = E 0 
A G 


24,132], Define <f>(z) = S(z) - A(E(z) 


E 0 ) for some 

which plays the role of a Lagrange multiplier. The states z^ that 


maximize the entropy under the constraint E(z) = E$ are solutions to 


^ = 0, 

E(z) = E„, 


which is equivalent to 



Eii 


(13) 
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If there are more constraints then the functional and Eq. (1T5D need to 
be modified accordingly. Furthermore, since is also a non-decreasing 

function of time, it follows that all solutions z t of the GENERIC equation 
fjl2]) converge to as £ —>■ oo. We refer to [iH, Property 3] for more detailed 
discussions. 


2.2. Generalised GENERIC 

In 32], the author introduced a generalization of the GENERIC frame¬ 
work to the case where the dissipative part of the GENERIC evolution (I12p 
does not depend linearly on gradient of the entropy ||. To describe the 
setting in [32[, one needs a dissipation potential /C such that for all z, the 
function /C(z; •) is convex and nonnegative and satisfies /C(z; 0) = 0. The 
convex sub differential of /C(z; •) is defined by 

^/C(z,0 = {^|/C(z;0>/C(z;0 + (f-e,^) for all £}. 

If £ H* /C(z, £) is differentiable, then the sub differential c^/C(z,£) is the same 
as the derivative of /C(z, £). 

The generalized GENERIC is then given as the differential inclusion 


9 i z6L(z)5’+^(z; iS(z) 


6z ? V ’ 5z 

and the degeneracy condition (fill) is replaced by 


Vzez, £gZ*, Ag 


L(z)^^ = 0, and /c(z;£ + 


5 z 


5z 


(14) 


= K(z;{). 


In particular, the energy is preserved and the entropy is non-decreasing along 
solutions of the generalised GENERIC evolution (1T4l) . We refer to [323, Sec¬ 
tion 2.5] for more detailed discussions. 

2.3. GENERIC formulation of the heat equation and the kinetic Fokker- 
Planck equation 

The heat equation (P) can be written as [l] 


where 


d t p= M(p)^, 
dp 


M(p)f = —z/div(pV£), S(p) = - / plogp, 


(15) 

(16) 
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which is a special instance of the GENERIC framework, in which the con¬ 
servative part is absent. Note that this GENERIC structure coincides with 
the Wasserstein gradient flow formulation of the heat equation BH0. 

To write the kinetic Fokker-Planck equation (J2]) in the GENERIC frame¬ 
work, in 18], the authors introduced an axillary variable e, depending only 
on t, so that the total energy e(t) + f (j- + V(q))p t (q, p)dqdp is conserved, 
leading to 

a r „2 

(17) 


d , 

77 e = 7 / 
dt JM? d 




m 


The coupled system (l2l)- (fT7l) was then shown to be a GENERIC with the 
building blocks 


2 ,p{^ 2d ) X 


Z = v : 

z = (p,e), 

L = L (p, e) = f q P Q 


E(p, e) = 


p 

2m 


+ V{q) )p(q,p)dqdp + e, 


S(p, e) = -6 I p(q, p) log p(q, p)dqdp + e, (18) 

M, 




where the operators defining L and M are given, upon applying them to a 
vector (£, r) at (p, e), by 


L PP £ = divpJV£, 

Mpp^ = - divp pV p £, 

= - / — • p(dqdp), 

Jr™ m 


1 pe r = rdiv p (p£), 


M PP r = r 


/ JL 

/jj 2 d m 2 


p(dqdp). 


Here 


J = 


I 0 


(19) 


is a 2d- dimensional anti-symmetric matrix, and Vo iP (M 2a! ) is the space of 
probability measures on with bounded second p-moments: 


V, 


2, p( 


t>2d\ 


) := \p G 'P(M 2d ) : / p z p(dpdq) < oo 





3. GENERIC formulation of the relativistic heat equation and of 
the relativistic kinetic Fokker-Planck equation 

In this section, we cast the relativistic heat equation (J3]) and the rela¬ 
tivistic kinetic Fokker-Planck equation (J5]) into the GENERIC framework 
and compute the stationary solution of the latter. 


3.1. GENERIC formulation of the relativistic heat equation 
We define the dissipation potential /C by 


£(p;0 = v / pip*(X7£)dx, 


where 


<P K z ) = ~2 


1 + 


v“ 


-1 


hence V<p*(;?) = 


1 + 




It follows that 

<9g/C(p, £) = — z/div(pV<p*(V£)) = —I'div | p 

'i + Slv£| 2 , 

Hence the relativistic heat equation (J3]) can be written as 

Vlogp 


dtp = v div p 


= ^/c(p, 


1 + ^j|V logp| 2 
where the entropy S(p) is the negative Boltzmann entropy 

S(p) = - P logp. 


ss 

Sp 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


It is straightforward to check that /C(p, •) is convex, nonnegative and satisfies 
/C(p, 0) = 0. Hence, Eq. (1221) with the building blocks (1231) is an instance of 
the generalised GENERIC evolution in which the conservative part is absent. 
Note that this generalised GENERIC structure coincides with the gradient 
flow structure of the relativistic heat equation shown in 31]. In addition, 
since lim^oo p*(z) = ^j-, at least formally one obtains the GENERIC struc¬ 
ture (|T5j) - (rr6]) of the classical heat equation as c —> oo. We comment on this 
more in Section [2 
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3.2. GENERIC formulation of the relativistic kinetic Fokker-Planck equation 
In this section, we show that both variants of the relativistic kinetic 
Fokker-Planck equation, i.e. Eq. (J3J) and Eq. (15]h can be formulated in 
the GENERIC framework. We define the relativistic Hamiltonian H by 

H(q,p) := c\/m 2 c 2 + |p | 2 + V(q), (24) 

and the diffusion matrix by 


D := 



_ me 

m 2 c 2 +\p\ 


rA J + 


p®p \ 
m 2 c 2 / 5 


for Eq. (SD, 
for Eq. (|5]). 


(25) 


For each p G M d , D(p) is semi-positive definite since for all £ G it holds 
that 


D(p)£-Z = 


, 2 2 7 I I2\ y](™ 2 C%+m)66 

m[m z c z + \p\ z ) 

v 11/ ij = 1 


1 


c ||£|| + (p ■ £) ) > 0. 


m{m?c 2 + |p| 2 ) 

We compute V p i7 and D(p) ■ V P H as follows. We have 

cp 


and 


D{p) ■ V P H = 


V P H = 


me 


\Jm 2 c 2 + |p| 2 ’ 


/ + 


p®p 


cp 


yjm^c 2 + |p | 2 V m 2 c 2 J yjm 2 c 2 + |p | 2 

= / 22,1 m (ra 2 c 2 / + P ® P) ' P 

m[m z c z + |pp) 

= 7 2 2,1 i 2 . (^ 2 C 2 + |p| 2 ) P = —• 

m(m z c z + |pp) m 

Hence both Eq. (J4]) and Eq. ([5]) can be re-written as follows 


d t p = div (JViPp) + 7 divp [O(0V p p + pV p H)\ , (26) 


with the corresponding matrix D, and the matrix J is given by 


J = 




5 
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being a 2d- dimensional anti-symmetric matrix. Let p t be a solution of (125]) . 

Suppose that V is non-negative and that the initial data p 0 satisfy 
,/ E 2 d H(qpp) p 0 (dqdp) < oo. We now show that f R2d H(q,p) p t (dqdp) is finite 
for any t > 0. Indeed, let us compute 


d 

dt 


H(q,p) pt(dqdp) 


H(q,p)d t p t (dqdp ) 


H(q,p) ( di v(JVHp t ) + 7 div p [B(9X7 P p t + p t S7 v H)} 


(*) 


= -7 / ©VpiL • [/yVpiL + 0 V P pt] dqdp 


= —7 / [DVpiL • X7 p Hp t + 0BV p ii • V p pt] dqdp 


= -7 / ©Vpif ■ V P Hp t + 7 9 diVp(DVpiL) p 4 , 


(27) 


where in (*) we have used the fact that J is anti-symmetric. Note that 


BV p ii ■ V p i7 


IV i /| 2 = c2p2 _ 

| V p | m 2 c 2_|_„2 5 

o 

cpj 


\yj m 2 c 2 +p 2 


for Eq. (HD 
for Eq. ©. 


and 


diVp(DVpii) 


A V H = 

diVp(^) 

V v m' 


d 

ra 2 c 2 +p 2 

d 

m 5 


(m 2 c 2 +p 2 )2 


T I < 


m ’ 


for Eq. (HD, 
for Eq. (ED- 


Substitute this computation to (I 2 ?D we get 


d 

dt 


H(q,p ) pt(dqdp) = -7 [ BVpii ■ Vpiipt + 70 


divp (BVpii) p t 



2 2 

m 2 c 2 P + p2 Pt + lO f R2d A pH p t , 

cp 2 „ 1 

-A + W 


rtyt 


m 2 c 2 +p 2 


for Eq. (HD 
for Eq. (ED- 


( 28 ) 
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Since ID'VpH ■ V P H > 0, the above computation implies that 


H(q,p) p t (dqdp) <—. 

(it /na>2d 777/ 


( 29 ) 


Applying Gronwall’s inequality, we obtain 
[ H(q,p) p t (dqdp) < —T + f 

J K 2d m J R ; 

for all t > 0 and some C > 0. 

As a consequence, since V > 0 , we have 


H(q-p) po(dqdp) < C, (30) 


c 2 p 2 


m 2 c 2 + p- 


<c 2 , 


cp 


( m 2 c 2 + p 2 ) 2 


C\Jm 2 c 2 + p 2 H(q,p) 


4 m 2 c 2 


Am 2 c 2 


cp 


m\Jm 2 c 2 + p 2 m 


< —a Jm 2 c 2 + p 2 < — H(q,p). 


m 


(31) 

(32) 


It also follows from the above computation that f R2d H(q,p)p t (q,p) dqdp is 
not preserved. Similarly as in Section [2731 we introduce an excess variable e, 
depending only on t, such that 

TT^) + 37 / H(q,p)p t (q,p) dqdp = 0 , 
dt dt J-^2d 

which, from (l27lh results in 

d 
dt 


e(t) = 7 / DVpi/ • V p Hp t - 7 9 div p (OV p i7) p t 


+ 


'I'JrM m 2 C 2 +p 2 A + C 3^/ R 2d ( yV2C 2 +P 2 ( m 2 c 2 +p 2)5 

7/lE 


cp" „ _ Tfrf 

m 2 c 2 +p 2 ^ m ’ 


) p t , for Eq. 
for Eq. 
(33) 


Note that the right-hand side of (j33j) is well-defined due to (l30]h (13D and 


We emphasize that the coupling is only one-direction: Eq. (1331) is slaved 
to Eq. (1261) . In other words, the addition of the excess variable e is simply 
a mathematical technique, and does not change the original system. We 
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now place the coupled system 0260 - 03.31) in the GENERIC framework. The 
building blocks are defined as follows. 


Z = ^ 


d 2 d\ 


X 


E (p,e)= H(q,p)p(q,p)dqdp + e, 


z = (P,e), 

L = L(p, e) = f ™ Q 


S(p,e) = -9 I p(q,p)logp(q,p)dqdp + e, 

'U nn M, 


M = M(p, e) = 7 


<pp lvl pe 

M ep M ee 


(34) 

where the operators defining L and M are given, upon applying them to a 
vector (£, r) at (p, e), by 


Lpp^ = div (pJV£), 

Mpp^ = -diVp(DVp^p), M pe r = r div p (BV p H p), 

M ep £ = — / ©Vphf • X7 p £ p(dqdp)i M ee r = r / D V P H ■ V P H p(dqdp). 
J R 2d J R 2d 


We now verify that the above building blocks indeed give rise to a GENERIC 
formulation of the coupled system (I26j) - (l33l) . First of all, we compute deriva¬ 
tives of E and S 


JE(z) 

Sz 


H(q,p) 

1 


jS(z) 

5z 


- 0 (logp+ 1 ) 
1 


(35) 


Substituting these to definition of the Poisson operator and the dissipative 
operator, we obtain 


L( z )4E^) = (div(pJV-ff)j f 

It follows that the equation 

d t z = 


Mfz 


<SS(z) 


L(z 


<*E(z) 

“ 


= n 

M(z) 


diVp(DVpp) + diVp(DVphfp) 

f R 2 d DVpih • [ 9V p p + V p H]dqdp 


jS(z) 

5z 


is indeed the same as the coupled system f[26]) -fl33 |) . Next, we check the 
conditions put on the Poisson operator L. Let (£i> r i) and (£ 2 ,r 2 ) be two 
arbitrary vectors at (p, e). We have 


((£i,n), L(p,e)(£ 2 ,r 2 )) = ( 6 , L P p(p)£ 2 ) = f divpJV £ 2 

J R 2d 


V6•JV6p 


• JV6P = -((6,r 2 ),L(p,e)(6,r 1 )), 
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where we have used that J is anti-symmetric. The above equality shows that 
L is antisymmetric. The verification of the Jacobi identity (jBJ) is a lengthy 
but elementary calculation, which hinges on the fact that J is constant and 
antisymmetric, and hence is omitted. Now we verify the symmetry and 
positive-definiteness of M(z). Let (£1, ri) and (£2, r2) be two arbitrary vectors 
at (p, e), then we have 


((£i,c), M(p, e)(£ 2 ,r 2 )) 


= 7 


= 7 


= 7 


£1 div p (DVp£ 2 p) +r 2 / £1 div p (D V P H p) - n / D V P H • V p £ 2 p 


rir 2 / DVpiL • VpR p 


DV P £ 2 ■ V p £ 1 p-r 2 j BX7 p H • V p £i p + rq / £ 2 div p (DV p iL p) 


r 2 ri / DVpLJ ■ VpR p 


- / £ 2 diVp(DVp£ip)+ri / £ 2 div p (DV p iLp) - r 2 / DV p iL • V p £i p 


+ r 2 ri J DVpiL • V p i7 p 
((£2,^2), M(p, e)(£i, ri)), 


i.e., M is symmetric. It also follows from the above computation, by taking 
(£i,n) = (£2,^2) = (£,r), that 


((£, r), M(p, e)(£, r)) = 7 J ®(V P £ - rV p H) • (V p £ - rV p iJ) p > 0, 
i.e., M is non-negative. 

Note that the Poisson operator is the same as in Section 12.31 The coinci¬ 
dence of the L-operator was also recognized between the relativistic and 11011- 
relativistic imperfect fluids equations [34|. In addition, at least formally we 
have linic^oo Vpi/ = —, lim ( .^ 00 D(p) = /; hence we recover the GENERIC 
structure of the kinetic Fokker-Planck equation presented in Section 12.31 as 
c —y 00 . We also comment on this more in Section [4l 

Remark 1. The crucial property, which is true for both Eq. (|4j and Eq. 
P and makes the calculations in this section work, is that the drift term in 
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p-direction (i.e., the friction term) is exactly 7D VH. Note that this relation 
also holds true in the classical kinetic Fokker-Planck equation and is known 
as the fluctuation-dissipation (or Einstein) relation. It is this relation that 
allows us to write all these systems into a common formulation as in (1261) . 
We note that this relation is more general than the one mentioned in [12L 
Eq. (13)] that concerns only the coefficient of the drift term. 


3.3. Stationary distribution of the relativistic kinetic Fokker Planck equation 
Next we seek stationary solutions of the coupled system (l26l) - fl33l) using 
m- A stationary solution = (p^, e^) of (f25]) - (153|) maximizes the entropy 
S(z) under the constraints that E(z) is a constant and that J R2d p(dqdp) = 1. 
Define <h(z) = S(z) — Ai(E(z) — E 0 ) — A 2 (/ R2 d p(dqdp) — 1) for some Ai, A 2 G M 
then (poo,eoo) satisfies the following equation 

f¥ = o, 

{ E(z) = Eo, 

{ J R 2 d p(dqdp) = 1. 

From (l35]h this system of equations is equivalent to 


'-0(log p{q,p) + 1 ) = Ai H(q,p) + A 2 , 

1 = Ai, 

J R 2dH(q,p)p(q,p)dqdp + e = E 0 , 
J r m p(dqdp) = 1 . 


Solving this system we obtain p 00 (dg,dp) = Z 1 exp ( — 1 H(q,p)) dqdp = 


Z 1 exp (-l(<V?n 2 c 2 + \p\ 2 +V(q ))) dqdp and = E 0 -/ R2d H(q,p)p 00 (dq, dp), 

where Z is the normalizing constant 


Z = 


exp 


( - i(c\/m 2 c 2 + |p| 2 + V(q))j dqdp. 


Note that p 00 is exactly the relativistic Maxwellian distribution. In other 
words, we obtain the relativistic Maxwellian distribution as the stationary 
solution of the relativistic kinetic Fokker-Planck equation as a consequence of 
the GENERIC formulation. This stationary solution has also been obtained 
in 0, 0 0 using different arguments. Our approach works simultaneously 
for both Eq. (@J and Eq. ([5]) and interprets why they have the same station¬ 
ary distribution. 
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4. Conclusion and outlook 


We have unified the relativistic heat equation and the relativistic kinetic 
Fokker-Planck equation into the GENERIC framework analogously to their 
non-relativistic counterparts. We also have established that the Maxwellian 
distribution is the unique stationary solution of the relativistic kinetic Fokker- 
Planck equation. A possibility for further research would be to study asymp¬ 
totic limits of the relativistic kinetic Fokker-Planck equation such as massless 
limit or the Newtonian limit, i.e., when the speed of light goes to infinity, 
using the GENERIC structure. It has become clear, see e.g., [40, Hi], [o] an 
a recent paper 33], that the formulation of dissipative equations via the en¬ 


tropy functional and the dissipation potential is compatible well with asymp¬ 
totic limits. This is because many techniques in calculus of variations, such 
as Gamma convergence, can be exploited. A similar framework for fully 
GENERIC systems has been recently proposed in 18]. We expect that these 


asymptotic behaviour of the relativistic equations can be obtained using the 
framework in the mentioned papers. We will return this issue in future pub¬ 
lications. 
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